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Abstract 

We show that the detection of neutrinos from a typical gamma ray burst requires 
a kilometer-scale detector. We argue that large bursts should be visible with the 
neutrino telescopes under construction. We emphasize the 3 techniques by which neu- 
trino telescopes can perform this search: by triggering on i) bursts of muons from 
muon neutrinos, ii) muons from air cascades initiated by high energy gamma rays and 
iii) showers made by relatively low energy (~ 100 MeV) electron neutrinos. Timing of 
neutrino-photon coincidences may yield a measurement of the neutrino mass to order 
10~ 5 eV, an interesting range in light of the solar neutrino anomaly. 



1. Introduction 



The origin of gamma ray bursts (GRBs) is arguably astronomy's most outstanding puzzle [T]]. 
Contributing to its mystery is the failure to observe counterparts in any other wavelength 
of light. It should therefore be a high priority to establish whether GRBs emit most of their 
energy in neutrinos M, [5], |J as expected in the (presently favored) cosmological models. 

It is not the purpose of this paper to study the modelling of GRBs. We will consider 
two cosmological scenarios: ultra-relativistic fireballs [||] and cosmic strings || and reduce 
their predictions to dimensional analysis, omitting details which represent at best unfounded 
speculations. After imposing experimental constraints on the dimensional analysis, it suffices 
to quantitatively frame the question of neutrino emission. The "experimental facts", which 
will later constrain our model parameters, can be encapsulated as follows^: i) there are 
about 100 bursts per year with an average fluency in photons of F 1 > 10 -9 Jm" 2 , ii) they 
are concentrated, on average, at a redshift of z ~ 1, iii) some bursts last less than 10 s, and 
iv) they do not repeat on a time-scale of 1 year or less. Our predictions will be presented in a 
form in which they can be scaled to fit varying interpretations of the experimental situation. 
Our interpretation of the observational situation, as well as the models presented, seem to 
be currently favored, although there are some dissenters. For example, some advocate that 
the origin of GRBs can be traced to an extended halo population of neutron stars. However, 
the predictions of such models for neutrino emission may in the end differ only slightly, since 
the reduced luminosity, compared to large-redshift sources, is compensated for by a reduced 
distance to the source. 

Our results can be summarized as follows. The detection of typical GRBs requires 
kilometer-scale neutrino telescopes. GRBs provide us with yet another example of Nature's 
conspiracy to require kilometer-size detectors for exploring our science goals 0, from dark 
matter searches to the study of active galaxies. Rare, large bursts may however be within 
reach of the present experiments. Our results will demonstrate that non-observation will 
lead to meaningful constraints on the models. In particular, it is unlikely that cosmic string 
models can escape the scrutiny of the detectors presently under construction, because they 
predict a fluency in neutrinos which exceeds that for photons by a factor of order 10 8 or more. 
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Furthermore, we will emphasize the 3 techniques by which neutrino telescopes can search 
for GRBs. All detectors such as the DUMAND and NESTOR deep ocean experiments, 
can search for short bursts of high energy muons of i/^-origin. Sensitivity is good, i.e. atmo- 
spheric backgrounds small, because the signal integrates over very short times and does not 
have to be searched for; one looks at times given by the gamma ray observations. The shal- 
lower detectors like AMANDA and Baikal can also search for the muons made in air showers 
initiated by TeV gamma rays|| of GRB origin. Finally, AMANDA can use its supernova 
trigger |[10| to identify excess counting rates in the optical modules associated with a flux of 
MeV-GeV u^s for the duration of a gamma ray burst. 

It has not escaped our attention that the observation of coincident bursts of neutrinos 
and gamma rays can be used to make a measurement of the neutrino mass. The mass is 



determined from the time delay td by simple relativistic kinematics with m v = E u ^j2ctd/ D. 
With td possibly of order milliseconds, distances D of thousands of Megaparsecs and energies 
E v similar to that of a supernova, neutrino observations from GRBs could improve the well- 
advertised limit obtained from supernova SN1987A by a factor 10 6 . The sensitivity of order 
1CT 5 eV is in the range implied by the solar neutrino anomaly. The measurement would be 
greatly facilitated by the fact that, unlike for rare supernova events, repeated observations 
are possible. 



2. Accelerator I: The Relativistic Fireball Scenarios 

Although the details can be complex, the overall idea of fireball models is that a large 
amount of energy is released in a compact region of radius R ~ 10 2 km~ cAt. The shortest 
time-scales, with At of order milliseconds, determine the size of the initial fireball ||. Only 
neutrinos escape because the fireball is opaque to photons. In GRBs a significant fraction of 
the photons is indeed above pair production threshold and produce electrons. It is straight- 
forward to show that the optical depth of the fireball is of order 10 13 ||. It is then theorized 
that a relativistic shock, with 7 ~ 10 2 or more, expands into the interstellar medium and 
photons escape only when the optical depth of the shock has been sufficiently reduced. The 
properties of the relativistic shock are a matter of speculation. They fortunately do not 
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affect the predictions for neutrino emission. 

For a ffuency F = 10~ 9 J m~ 2 and a distance z = 1 the energy required is 

E y — 2x 10 51 erg ( ) ( -) , (1) 

7 S V 4000 M P C / U0- 9 Jm-V V; 

using E 1 = 47rD 2 F. The temperature T 7 is obtained from the energy density 

p = El = -haT A , (2) 
H V 2 K J 

where h represents the degrees of freedom (/i 7 = 2 and h u = 2 ■ 3 ■ | for 3 species of neutrinos 

and antineutrinos) , V the volume corresponding to radius R and a = 7.6 x 1CT 16 J m -3 K~ 4 . 

We find that 

U1 / E 7 \ 1/4 A00km\ 3/4 

T 7 = 8MeV - — — y - — . 3 

7 \2 x 10 51 ergy \ R J w 

For neutrinos 

For a merger of n-stars, for instance, the release of a solar mass of energy of 2 x 10 53 erg 
implies a total energy emitted in neutrinos ~ 10 2 E 1 . The 7's are most likely produced by 
bremsstrahlung of electrons from vv annihilation. The actual predictions for the energy and 
time structure of the photon signal depend on the details of the shock which carries them 
outside the opaque fireball region of size R. The data suggest that the structure of these 
shocks is complex. Neutrinos, on the contrary, promptly escape and carry direct information 
on the original explosion. From we obtain T v ~ 2.5T 7 ~ 20 MeV. Using this and a 

total neutrino energy in the fireball of lO 2 ^ we obtain 

At bs = 0.3 msec ( — ) . (6) 



VlOO km 

The neutrino fluency is obtained from E utot / '(4ttD 2 



E utot \ /65 MeV\ /4000 Mpc^ 2 



iv; = io 4 Hi-M — =^ — r (7) 

^2 x 10 53 erg/ V E v ) \ D ) K J 

or more than 10 57 z/'s at the source. Notice that this prediction is rather model-independent 
because it just relies on the fact that a solar mass of energy is released in a volume of 
100 kilometer radius which is determined by the observed duration of the bursts. 



Although the ~ 100 MeV-neutrinos are below the muon threshold of high energy neutrino 
telescopes, the v e will initiate electromagnetic showers by the reaction [p e + p — > n + e + ) 
which will be counted by the AMANDA supernova trigger. 

A supernova with properties similar to those of SN1987A can cause a 10 second burst of 
neutrinos in the AMANDA detector with E u ~ 40 MeV. They produce positrons with, on 



average, half that energy. Detailed simulations [10] of the supernova signal in the AMANDA 



detector have shown that each photomultiplier tube (PMT) has a seeing radius d ~ 7.5 m 
for 20 MeV positrons. The number of events per PMT is given by 

#N, obs ~ N^d 2 ) . (8) 

The last factor estimates the probability that the v e produces a positron within view of the 
PMT. Here 

Kl = ^Apa E 2 u (9) 

with 

a = 7.5 x 10~ 40 m 2 MeV" 2 . (10) 

A is Avogadro's number and p the density of the detector medium. One should not forget here 
that the dependence of the cross section on neutrino energy is linear rather than quadratic 
above ~ 100 MeV. 

We have checked by Monte Carlo [|TTJ that the seeing volume scales linearly in the energy 
of the positron, or neutrino, up to TeV energies. Eventually the radius will cease to grow 
due to attenuation of the light. With absorption lengths of several hundred meters [|12j this 



upper limit is outside the range of where we will apply (§). Therefore, the event rate for 
GRBs is given by (|8]) with d = 7.5 m ( 20 MeV ) ^ ' Here 30 MeV is the positron energy which 
is, on average, half the neutrino energy given by Eq. (^). 

Can this signal be detected by simple PMT counting? Signal S, noise N and S/ VN, for 
an average burst, are given by 

S = 10- 3 events ( ( ( ^) > (") 



,5 x 10-V V20 cm) V 200 J 

N = 60 events (-**-) , (12) 

V0.3msec/ VlkHz7 V 200 / v ; 

S/VN = l0 ^(^Y 1/2 (Pl^) 2 (13) 

1 VI kHz/ V20cm/ V 200 / V ; 



AMANDA has been chosen for reference with 200 PMTs with a diameter -Dpmt of 20 cm 
and a background counting rate of roughly 1 kHz. With such low rates in millisecond times, 
observation obviously requires a dedicated trigger. 

Obviously the event rate for an average burst is predicted to be low. We will argue 
nevertheless that observation is possible and clearly guaranteed for kilometer-scale detector 
with several thousand PMTs. First, the parameters entering the calculation are uncertain. 
The event rate increases with neutrino energy as because of the increase of the PMT 
seeing distance d and the neutrino interaction cross section With increased energy the 
average burst may become observable. Individual burst can yield orders of magnitude higher 
neutrino rates because of intrinsically higher luminosity and/or smaller than average distance 
to earth. For example, a burst 10 times closer than average and 10 times more energetic is 
observable with a significance of well over 10 a in the exisiting AMANDA detector. Given the 
uncertainties in the model and its parameters as well as the chaotic nature of the phenomenon 
(there is no such thing as an average GRB), this event represents a plausible possibility. 

As demonstrated by the 7-ray observations, the structure of the shock producing the 
gamma rays is complex. The interaction of multiple shocks can also produce neutrinos on 
other time-scales and with different, sometimes much higher, energies^. So one should 
have an open mind when searching for bursts. This is underscored by the rather different 
predictions obtained from string-type models, which we discuss next. 

3. Accelerator II: Cosmic String- Type Scenarios 

The dimensional analysis relevant to accelerators such as cosmic strings is synchrotron emis- 
sion from a beam of ultra-relativistic particles. The time of emission is now given by 

Ail ab = — . (14) 

Here L is the size of the accelerator and 7 = /saturation /-^ is a ratio of electric currents, which 
is some large number. One main difference with the previous scenario is that the emission 
is relativistically beamed in a solid angle of size 7 -2 . The idea is that when accelerated 
currents reach a value /saturation it is energetically more favorable to radiate away the mass 
of the accelerating cosmic source, rather than sustain the high current. This happens for 
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instance at cusps in oscillating loops where the current becomes, theoretically, infinitely 
large. A mass \i per unit length L is radiated away in a time At. In dimensionless units, /i 
is, 

e = /4- (15) 
c z 

A dimensional estimate for L, the size of the cosmological accelerator, can be made as follows. 
The time over which a cosmic accelerator loses mass is clearly proportional to Lj\i or, in 
correct units, L/ec. We equate this to the only time in the problem: the lifetime of the 
universe at the redshift of the accelerator, 

— = £ 7: \ 3/2 , (16) 
ec (1 + z) i i l 

where cto = 6 x 10 27 cm and the proportionality factor £ = 1. So L = £ecio/(l + -z) 3 ^ 2 and 
we can now calculate the duration of the burst 

At obscrv = (1 + z)At Uh = (1 + z )Jl = io 17 £ 4t seconds . (17) 

In the accelerator frame (comoving frame) 

(e/KT 11 ) 

At com S £ ^ f seconds , (18) 

( 7 /10 3 ) 2 V ; 

The choice of units will become clear further on. The energy loss per unit length is indepen- 
dent of e with 

(19) 



^ = 48xl0 33 m'jm-s- 1 



At com £ VlOV 
A fraction ?7 7 is radiated away in 7-rays. 

The above equations are valid for cosmological strings or loops of false vacuum in grand 
unified theories. Near cusps in oscillating loops the particle currents become very large, 
creating a situation where the energy density exceeds that of the topological defect and 
the energy is released in a short localized burst of radiation. In string models there is a 
proportionality factor multiplying the r.h.s. of fll6|) which is of order £ = 10 3 rather than 
unity; see e.g. Ref. 0. From now on we will include this factor, so that our results can be 
directly compared to these models. 
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Imposing the "experimental facts" , listed in the introduction, on the dimensional analysis 
(with £ = 10 3 ) yields the following constraints |6|: 

10 2 < 7 < 10 5 
1(T 12 < e < KT 11 (20) 
1(T 10 < r? 7 < KT 9 

The critical result here is that to accommodate the time-scales as well as the fluencies in 
a large redshift source of this type, the fraction of energy loss into gamma rays is actually 
very small, 10~ 10 to 10~ 9 . Theoretical argument s|| lead to the expectation that most of the 
energy is radiated into z/'s. This fits well with the observational fact that the missing energy 
is not emitted in any other wavelength of light. 

Before proceeding it is important to point out that the small fraction of the burst energy 
going into gamma rays is not a surprise. Cosmic strings belong to the class of highly inef- 
ficient models in which the whole accelerator is boosted by a Lorentz factor 7. In contrast, 
conventional fireball models describe a collisionless shock of protons which carries kinetic 
energy far outside the opaque fireball where it is transformed into a burst of photons. 

A fraction 77" 1 is radiated into z/'s of energy E uo ^, s . The flux for a typical burst is 

^1^! (21) 

'/7 '-'v obs 

or 

N. per crf = 10 * (^y\ T ^y\ w I f -;) (22) 

during a time 

^ = 1 sec 3 



At ohs = - At com = 1 sec ( 77— n~ ) • (23) 



Here 



E uohs = ^3.15T ucom . (24) 

The thermal emission of the neutrinos in the accelerator frame follows a Fermi-Dirac distri- 
bution with temperature T„ com . We will estimate it next following Ref. [Q. 

Consider an accelerator segment of loop of length L and radius R. Assume black body 
radiation off its surface and apply the Stefan-Boltzmann law in a comoving frame. Using (|T9|) . 



/iC 
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L = (2tt RL)(aT* com ) , (25) 



where a is the Stefan-Boltzmann constant. We obtain 



For a cosmic string R = /saturation /Her, where H CT the critical field strength, /saturation was 



calculated by Witten|13|, and is typically 

i(T 8 m< R < l(r 6 m. (27) 

The possibilities covered by this class of models range from thermal supernova-type en- 
ergies to TeV-neutrinos. For illustration, we show results for a low and high energy neutrino 
scenario. 

7 = 10 2 R = KT 6 E uohs = 560 MeV 
2 x 10 7 < N uohs < 2 x 10 8 per cm 2 
10 2 < At obs < 10 3 seconds 

or 

7 = 10 5 R = 1(T 8 E vohs = 60 TeV 
2 x 10 2 < N uohs < 2 x 10 3 per cm 2 
0.1 < At bs < 1 /xsec 

Suppose neutrinos with E uohs ~ 40 MeV produce electrons in the detector with energy 
(1 - (y})E u , or about 20 MeV, just like SN1987A would have produced in AMANDA. We 
calculate a flux of 5 x 10 8 per cm 2 in a rather long burst. We know from the supernova 
analysis that each PMT has a seeing radius d ~ 7.5 m in this case. The number of events, 
given by (|8|), is 10 per PMT for a typical, average burst. This is 10 times smaller than a 
supernova, but the GRB data indicates that we have 100 shots per year and there should be 
some big ones. Models suggest searches over > 1 sec intervals, maybe up to 1000 sec. Also 
notice that event rates grow with energy as ad 3 /E. Both d, a grow with energy. The signals 
should be spectacular for Z^obs values of hundreds of MeV or more. 

An extreme example on the high energy end yields ~ 10 2 neutrinos of tens of TeV energy 
per cm 2 in periods -C 1 sec. In this scenario, the secondary muons can be detected and 
reconstructed. This allows one to both count the neutrinos and reconstruct their direction 
with degree-accuracy. The event rates are now given by|§: 

N ev ents = N„ Area P v ^ , (28) 



- "^="( l0 * 42m2 (S?)( 5m (^)' (29) 

Here Py-*^ is the probability that the neutrino interacts and spawns a muon that reaches the 
detector; it is proportional to the density p of the detector medium, the neutrino interaction 
cross section <j v and the muon range For E^ ~ \E V ~ 30 TeV and p = j^A per cm 3 we 
have P„^ M = 10~ 3 or 10 5 events for a detector as small as 100 m 2 area detector! 

Therefore, bursts associated with topological defects are unlikely to escape the scrutiny 
of both the supernova and the muon trigger. In part of the parameter space one should be 
able to rule out the cosmological models even for average bursts. In other regions, one can 
constrain the models only from a search for energetic bursts. 

4. Detecting 7-Rays with Neutrino Telescopes? 

What about seeing 7-rays? Shallow detectors like AMANDA and Baikal detect secondary 
muons produced by 7-showers in the atmosphere. For a vertical muon threshold of 180 GeV, 
AMANDA should be sensitive to TeV gamma rays. The number of photons is calculated 
from the fluency F 7 by 

NJ> E) = -^-, (30) 

where a is the spectral index (a = 1 for Fermi shocks). For a = 1 and a fluency per burst 
of 10 -9 J m -2 we find that F~ = 10~ 2 In" 1 ! ^ 7max ) per m 2 per burst. There is a rather weak 

\ -^7 min / 

logarithmic dependence on the maximum and minimum energy of the photons in the burst. 
Notice that the TeV flux, even if it exists, is too small to be detected by satellite experiments. 
The maximum energy of GRBs is therefore an open question. It has been speculated that 
they may be the sources of the highest energy cosmic rays which implies a very high energy 
accelerator indeed. 

The muon flux produced by above gamma ray flux can be computed following Halzen 
and StanevPI: 



5 F 7 1 f cos 9E jm£kX \ f E^/ cos 9 



0.53 



NJ>E U ) ~ 2xl0- 5 ^- 7 , —— In 7max M/ . (31) 

M ^' cos9 (E^/ cosd) a+l \ 10£ M j\ 0.04 j v ; 

Here E^ is the vertical threshold energy of the detector, e.g. 0.18 TeV for the AMANDA de- 
tector. 9 is the zenith angle at which the source is observed. This parametrization reproduces 
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the explicit Monte Carlo results. 

We predict 10~ 6 muons per m 2 for an average burst, which can therefore be detected 
in a km 2 telescope! The probability that a 1 TeV 7 contains a detectable muon is about 
10~ 4 . We assumed here a burst in the 1 MeV to 10 TeV range and cos8 = 1. All this 
requires, of course, that the GRB flux extends to TeV energies. We do not know whether 
any do because satellite experiments have no sensitivity in this energy range. There is no 
atmospheric \i background in a pixel in the sky containing the GRB on a 1 second time scale. 
Big bursts may be detectable in the 10 m 2 detectors presently under construction. 
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